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We study numerically classical and quantum dynamics 
of a piecewise parabolic area preserving map on a cylinder 
which emerges from the bounce map of elongated triangu- 
lar billiards. The classical map exhibits anomalous diffusion. 
Quantization of the same map results in a system with dy- 
namical localization and pure point spectrum. 

PACS numbers: 05.45. Pq, 05.45. Mt 



As the Schrodinger equation is linear, it cannot posses 
exponential instability in time as do the classical equa- 
tions of motion. Quantum time evolution is dynamically 
stable 01 and all Lyapunov exponents are strictly zero. 
This does not contradict the correspondence principle as 
the quantum evolution still mimics classical exponential 
instability up to so-called log-time logl/?i. For larger 
times, the classical evolution develops details on scales 
smaller than h and the strict correspondence is lost. In 
some quantum systems there is another time scale, con- 
nected with the phenomena of dynamical localization. 
This is purely quantum interference effect and results in 
a quantum suppression of classical diffusion. For times 
smaller than this so-called break time t < t^, quantum 
evolution follows the classical diffusion. Afterwards, the 
quantum system "notices" it's discrete spectrum and the 
diffusion stops, resulting in a localized state. Since the 
break time is usually much larger than the log-time, 
exponential instability may not be relevant for dynamical 
localization, although all known examples of dynamical 
localization, from kicked rotors and maps |^,^ to quan- 
tum billiards Q, take place in classically chaotic systems. 

There are several reasons for studying systems with- 
out exponential instability. Most important is to gain 
insight into the old question, namely what are the nec- 
essary and sufficient properties of a system in order to 
display certain statistical properties. This question can 
be addressed in classical as well as in quantum mechanics. 
In this Letter we study a simple dynamical model with 
purely parabolic, non-chaotic and non-integrable, classi- 
cal time evolution. By means of numerical experiments 
and heuristic arguments we show that the model exhibits 
anomalous classical diffusion whose mechanism is essen- 
tially different from the one that is commonly observed in 
softly chaotic, i.e. KAM, systems. Furthermore, we show 
that the quantization of the same model displays dynam- 
ical localization, supporting a provocative claim that dy- 
namical localization can exist without classical chaos and 
without normal classical diffusion. 

Triangular billiards are among the simplest classical 



systems without exponential instability. It is worth men- 
tioning that a triangular billiard is equivalent to the mo- 
tion of three particles on a ring (provided all three angles 
are smaller than 7r/2) with their masses connected to the 
angles of a triangle pi . Here we consider a special asymp- 
totic case of an elongated triangular billiard with one an- 
gle 7 being very small where it was shown that in 
the limit 7^0 the bounce map, relating two successive 
collisions with the short sides of a triangle read 

Vr,+i = ^;„ + 2(u„ - [-u„] - /x(-l)[''"l) 

M„+i = M„ - 2i;„+i (mod 2), (1) 

where u, v are suitable position-momentum Birkhoff co- 
ordinates, ^ is a parameter measuring the asymmetry 
of two non-small angles, and is the nearest inte- 
ger to Un. Definition range for variables {u,v) is a 
cylinder. Furthermore, if compactified onto a torus 
(m (mod 2),w (mod 1)), the map (|^) can be rewritten in 
new variables x,y G [—1,1), with ?/„ = 2(— l)"(u„ + Vn ~ 
i) and Xn = (-l)"(u„ - 5), 
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where sgnx„ denotes the sign ±1 of Xn. Both maps, 
the map (|l]) and the so called triangle map (|^) are area 
preserving and (piece-wise) parabolic, therefore they do 
not exhibit exponential instability. We stress once again 
that the maps (|l]) and (^) are equivalent only if both 
are considered on the 2-torus. In the present paper we 
will focus on the classical and quantum properties of the 
map (|l|) on the cylinder. Let us first discuss its classical 
properties. 

If considered on the torus (^, three different cases of 
parameters must be distinguished: (i) If the parameter 
fi = k/2l and the initial ya = 2m/l are both rational (xq 
can be arbitrary irrational) then the dynamics can be 
reduced to a permutation on a finite integer lattice / x /, 
hence all orbits are periodic with period which divides P. 

(ii) The parameter ^ = k/2l is rational, but a;o and yo are 
irrational. In this case, the dynamics is pseudointegrable 
and the motion appears to be ergodic on the invariant 
curve which is a set of I helices {y — yo)l (mod 2) = 0. 

(iii) If the parameter /i as well as initial coordinates xq, yo 
are irrational then the motion appears to be ergodic on 
the whole phase space. However, the exploration of the 
phase space is very slow, namely it was found numerically 
that the number of different values of yt visited up to time 
t grows as ~ logi, i.e. weak ergodicity |^]. 

If we consider the same three cases in the original vari- 
ables u and V, that is the map (^ on the cylinder, the 
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case (i) results in a ballistic transport as generic periodic 
orbits of the triangle map now result in ballistic orbits. 
In pseudointegrable case (ii) one may focus on the motion 
along infinite invariant helices {{u+v)l (mod 1) — const). 
Numerical evidence suggests strongly that the motion 
along the helix exhibits anomalous diffusion (w„ — wo)^ ~ 
n", with exponent 1 < a < 2 which depends on pa- 
rameters ^,uo,vo. In weakly ergodic case (iii), we have 
'logarithmically slow' diffusion transverse to the helices 
due to weak ergodicity and 'faster than normal' diffu- 
sion along the helix. As a result we find numerically 
significant power-law anomalous diffusion in momentum 
variable (w„ — wq)^ ~ Therefore, the transport prop- 
erties of the triangle map (|l|) in cases (ii) and (iii) are 
quite similar, hence the irrationality of fi does not seem 
to play very important role, apart from determining the 
dimensionality of the invariant surface. From now on, we 
set in all numerical experiments ^ = {^{Vb— 1) — e~^)/8 
(same as in [^), but note that results are practically in- 
dependent of this particular value. 
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FIG. 1. Scaling function f{v/t°'^'^). Scaled density distri- 
butions for three different times t — 6000, t = 12000 and 
t = 18000, full, dashed and dotted curves, respectively, are 
shown and all fall on the same curve. Averaging over 20 • 10^ 
initial conditions is performed. Asymptotic solid line is expo- 
nential function. 

If we denote by p(v, t) the density of a swarm of trajec- 
tories along the momentum axis as a function of time t, 
J pdv — 1, then the dispersion a^{t) = J v'^p{v, t)dv (x t" 
grows faster than linear and defines exponent a (see 
the inset in fig. ||). For our map (|l|) the exponent 
of anomalous diffusion is approximately a — 1.81 ± 
0.01 and only weakly depends on p. In addition to 
that, as we show in figj^, the momentum-density start- 
ing from initial p{v,Q) = S{v) satisfies a scaling law 
p{v,t) = f~"/^/(u/t"/^) where the scaling function f{x) 
has asymptotically (|x| —^ oo) exponential tails. Anoma- 
lous diffusion is usually explained in terms of Levy walks 
which emerge as a consequence of some hierarchical struc- 
tures in phase space that exhibit scaling. For smooth 
Hamiltonian systems, these structures are stable KAM 



islands with sticky boundaries. Orbit wanders through- 
out the chaotic sea, but now and then it gets trapped 
around one of those islands, causing temporary ballistic 
motion. However, it seems obvious that successive ballis- 
tic fiights are uncorrelated as the orbit spends the mean- 
time in the chaotic sea with positive Lyapunov exponent, 
justifying a random walk approximation, which gives re- 
lation between the exponent of anomalous diffusion and 
the distribution of flight lengths for individual ballistic 
periods To our knowledge, all so far studied cases 
of anomalous diffusion involved systems with Lyapunov 
chaos, except recently studied "square-ring" billiard 
For systems without exponential instability the connec- 
tion between Levy flights and the exponent of anomalous 
diffusion cannot be so simple. This is clearly conflrmed 
by our experiments. The distribution of individual flight 
lengths has been found to fall as (Aw)^'^ with [3 w 2.9, 
whether according to the random walk approximation 
this exponent should be/3 = 4 — a~2.2. The reason for 
the discrepancy are strong correlations between consecu- 
tive flights. We have computed the correlation function 
between flight lengths (Aw(t)Aw(O)) and found that it 
falls only very slowly, perhaps as t^" with the exponent 
around k « 0.1. 

We note again the fact that the map (|^) exhibits 
anomalous diffusion also for a pseudointegrable case of 
rational p, and irrational j/o- Irrational j/o can be suc- 
cessively approximated by rational initial momenta j/o ~ 
2mi/li, e.g. by its continued fraction approximations, for 
which the orbit is ballistic. Since the system is parabolic, 
such approximate orbit will follow an exact orbit up to 
times of order t ^ U, which is typically smaller than the 
period of ballistic orbit compactified on a torus ~ tf. 
Hence there is no contradiction between the ballistic 
spreading for the rational approximates and the anoma- 
lous spreading for the irrational initial conditions. The 
question of explaining anomalous diffusion can perhaps 
be reduced to the number theoretic study of periodic or- 
bits. Note also that anomalous diffusion implies fractal 
singular continous spectrum of classical time evolution 

We now turn to the quantum version of the map 
(|l|). Quantization goes along the standard lines (e.g. 
Q), writing the Hamiltonian of the map (||) as H = 

T(v) + 5p{t)V{u), where 5p{t) = Y.n ^(^ ' kinetic 
part is T(v) = v^ = —d^/du^ and the potential is 

V{u) = i\u\+p-l/2f, for <1, (3) 

and is periodic in u with the period 2. Here we shifted 
the origin of u to the 1/2 with respect to (|l]) in order 
to have a symmetry V{u) = V{—u). Note that potential 
V{u) has a discontinuous derivative at integer values of 
u. In order to have a finite Hilbert space for numerical 
work, we compactify our cylindrical phase space to the 
torus {u (mod Q),v (mod P)), where Q = 2 and P is the 
periodicity (number of fundamental cells) in w-direction. 
Because we want to study diffusion, periodicity in the 
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momentum P will be set to some large integer, so that the 
phase space will effectively mimic a cylinder. Using peri- 
odic boundary conditions {u + Q\ip) = {u\^Ij), {v + P\iIj) = 
{v\ip), in position, momentum representation, respec- 
tively, Ti must satisfy QP/2ttTl = N, where N is an integer 
giving the dimension of the Hilbert space. There are nat- 
ural finite position and momentum bases satisfying the 
above boundary conditions, namely \uj) and \vk) respec- 
tively, with Uj = {j - N/2)Q/N,j = 0, 1, ... iV - 1 and 
Vk = {k-N/2)P/N, fc = 0, 1, . . . 7V-1. We perform quan- 
tum evolution li^it)) = U*|'0(O)) by decomposing the uni- 
tary propagator U = exp(— zT/?i) exp(— iV/?i) to diago- 
nal kinetic and potential part, with entries exp{—ivf /h) 
and exTp{—iV{uj)/h) respectively, where Fourier transfor- 
mation between the two bases {uj\vk) = N^^/'^e^^^'"''/^ is 
coded efficiently using FFT-algorithm. 
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FIG. 2. Anomalous diffusion in a classical map (|l|) (solid 
curve), and quantum cases for M — 100,50,25 (dot-dashed 
curves), and M = (lower inset). In the upper inset we 
show classical spreading (solid curve) and t" with best-fitting 
exponent a = 1.81 (dashed curve). 

Our main goal was to simulate time evolution and 
check quantum-classical correspondence with anomalous 
diffusion for short times and possibility of dynamical lo- 
calization for longer times. We argued earlier that lo- 
cal (in)stability of a classical system may not play a key 
role in dynamical localization. We know that in cases 
of normal classical diffusion one often finds pure point 
spectrum of quantum evolution, like e.g. in kicked rotor 
and that ballistic classical spreading, since it is typi- 
cally associated with regular structures in classical phase 
space, implies (at least partly) continuous quantum spec- 
trum. However, the case of anomalous classical diffusion 
is much less clear. At least there are no examples with an 
anomalous classical diffusion, that would exhibit dynam- 
ical localization, i.e. pure point quantum spectrum. It 
is therefore important to check the nature of a quantum 
spectrum for our map. We have done this along two in- 
dependent lines: direct time evolution in the momentum 
space and the solution of an eigenvalue problem. 




FIG. 3. Localized state after t = 10000 iterations averaged 
over At — 200 successive iterations, (l/At) X^J"*"^* |(w|V'(t))|^, 
is shown. Wave function (full curve) has a power-law cx 
tail (dashed curve). In the inset is in semi-log scale drawn 
only central, exponentially localized part of a wave function. 
Dashed curve is an exponential function. Everything is for 
P/N = 6765/349196 {M = 50). 

Following the time evolution of an initial momentum 
state {v\4>{0)) := S{v), we can check whether we have 
pure point spectrum or not. This has been performed 
using two different methods. The first method, which 
we call torus approximation (TA) has been described 
above and is a standard approach for dealing with quan- 
tum maps on cylinders. We have chosen P/N to be a 
continued fraction approximation of a noble irrational 
nh — 1/{M + (tgm), where ctqm — {V5+ l)/2 and M 
is an integer which determines the size of inverse Ti. By 
increasing M we approach the semiclassical limit while 
increasing P improves the accuracy of TA. However, we 
should note that within TA due to periodicity P in mo- 
mentum direction v we face significant aliasing of power- 
law tails in Fourier transform of a non-smooth potential 
V{u). Thus we have devised an alternative quantum evo- 
lution with an absorbing boundary (AB) in the momen- 
tum space. We truncate the Hilbert space to the set of 
N -\- 1 momenta Vk = ttM with k = -N/2, N/2, 
where now we can take any irrational value of h since 
we do not require periodicity in momentum space. The 
propagator for AB method is a finite nearly-unitary ma- 
trix Uk,i = 6{k)exp{—ivl/h))£^k-i where 9{k) 1 if 
\k\ < N/2 and 9{k) = otherwise. Now we calcu- 
late the Fourier components of the potential propaga- 
tor = (1/2) J^j^duexp{i7rku — iV{u)/h) exactly in 
terms of complex error function, that is with no alias- 
ing, however we may loose some quantum probability 
leaking from the system for large momenta. Control- 
ling the stability of the results, e.g. the quantum dis- 
persion a^{t) ~ {■i(j{t)\-v'^\ip{t)) and the stay probabihty 
ri{t) ~ {ip{t)\ip{t)) with increasing N we observe that AB 
method is typically much more efficient than TA, i.e. it 
converges for smaller dimension N . In general we found 
that results on a'^{t) of both methods agree up to some 
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time scale which monotonically increases with increasing 
A''. Beyond that time, o-^(i) of TA tends to overshoot 
due to ahased probabiUty at large momenta, while cr^(t) 
of AB remains stable. 

In fig.^ time dependence of a quantum dispersion 
cr^(t) is shown for several values of M together with 
the classical case. For cases M — 100, M = 50, 
and M = (inset) TA results are shown for P/N = 
6765/687446,4181/215815, and P/N = 121393/196418, 
respectively, which do not deviate significantly from the 
AB results within entire time range t < 10^. How- 
ever, for M = 25 AB result for N = 60000 satisfying 
T7(t) > 1 — 6 • 10~^ is shown instead, since we were not 
able to reach convergence of TA in the entire time range 
within accessible values of N. One can see that the quan- 
tum evolution follows the classical diffusion up to some 
time, after which dynamical localization sets in. The lo- 
calization is clearly demonstrated for small M, but for 
smaller Ti (e.g. M — 100) the localization length seems 
to increase very fast with M and the saturation could not 
be reached numerically within a reasonable time of com- 
putation. However, using standard heuristic arguments 
1^ , namely that the quantum spread roughly follows the 
classical ^ t" until the relaxation (break) time t = tB 
which is proportional to the density of 'operative eigen- 
states' ~ h~^a, we derive an asymptotic ?i-dependence 
of the localization length, a\{h) cx l/?i"/^^~"\ In our 
case, since a is quite close to ballistic value of 2 giving 
formally cri = oo, the localization length is predicted to 
grow very fast asymptotically, cti oc h~^'^ (x M^-^. Of 
course, we should be aware that the asymptotic scaling 
ai{h) is effectively not very accurate for not very small 
h and not very long times t where the classical spread 
increases with a smaller and varying effective exponent. 
In addition, localized states of time evolution {vk\ip{oo)) , 
as well as eigenstates of U, have an interesting struc- 
ture: an initial exponential decay with a rapid crossover 
to an asymptotic power law tail (~ v~'^) which is ex- 
pected [|l^ due to non-smooth C° potential (H) (see fig. 
g) . For smaller M the width of the central exponentially 
localized part of the wave function is smaller and the 
main contribution to a'^ comes from the spiky fluctuating 
power-law tails, as a results of which a'^{t) exhibits big 
fluctuations as a function of time (see insets in figs.||,H). 
On the other hand, in the semiclassical limit the central 
exponential part of the wave function becomes wider and 
consequently the tails become less important. 

The results of a stationary analysis of quasi-energy 
spectrum and eigenstates of the quantum propagator are 
fully compatible with the results of direct time evolution. 
We have diagonalized the matrix Ufc^; within TA and all 
eigenf unctions were found to be localized. We also ana- 
lyzed the spectral nearest neighbor spacing distribution 
P{s) which, as expected, is always found to have a Pois- 
son shape exp(— s) (see fig. ||). We must remark that by 
decreasing the size of the Hilbert space N at fixed h (or 
M), we obtain a transition to random matrix (Wigner- 



like) spectral distribution when P becomes comparable 
or smaller than localization length cr/. 
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FIG. 4. Level spacing distribution P{s) (full curve) for 
M = and P/N = 6765/10946. Dashed curve is a Pois- 
son distribution. In the inset we show a typical localized 
eigenstate of U for P/N = 28657/46368 (M = 0) peaked at 
vq — 10858.9 and the straight line gives the slope —4. 

In conclusion, we have demonstrated dynamical local- 
ization in a quantum system whose classical counterpart 
is a non-integrable, piecewise parabolic (locally stable) 
map and exhibits super-normal anomalous diffusion. The 
result implies that local dynamical instability of a clas- 
sical system is not necessary for dynamical localization. 
The mechanism of anomalous diffusion is here more com- 
plicated than in hierarchical chaotic systems, as the con- 
secutive Levy fiights are strongly correlated. 
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